Localized tendency for the superfluid and Mott insulator state 
in the array of dissipative cavities 
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The features of superfiuid-Mott insulator phase transition in the array of dissipative cavities is 
analyzed. Employing a kind of quasi-boson and a mean-filed approach, we show analytically how 
dissipation and decoherence influence the critical behaviors and the time evolution of the system. 
We find that there is a localized tendency, which could lead to the break of superfiuidity for a 
superfluid state and suppress the appearance of the long-range order form a Mott state. Eventually, 
a collection of mixture states localized on each site will arise. 
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One of the remarkable applications of coupled cavity 
arrays is to realize so called quantum simulators. Due to 
the controllability of atomic and optical systems, it could 
be useful to attack some unclear physics and to explore 
new phenomenon in quantum many-body systems [l|-|Jl. 
In particular, over the past years the experimental pro- 
gresses in controlling optical system and in fabricat- 
ing large-scale arrays of high-Q cavities j£, JLQ] make this 
potential application close to reality. However, the quan- 
tum optical system is typically driven by an external field 
and coupled to the environment [TTI . Il2j | , which bring the 
system out of equilibrium and profoundly affected the 
dynamics of interest (isl . . New important questions 
thus arise and need to be clarified, such as whether the 
link between the initial ideas of cavity arrays as quantum 
simulators and the realistically experimental conditions 
is hold, and how dissipation and decoherence would be- 
have in the system. 

In this paper, we propose answers to these questions 
by investigating the superfluid(SF)-Mott insulator phase 
transition in the array of dissipative cavities. Provided 
the external time dependence is much slower than the in- 
ternal frequencies of system, we show that there are still 
two fundamentally different quantum states, i.e. photons 
localized in each cavity(Mott-like) and delocalized cross 
the lattices(SF-like). However, a localized tendency does 
exist for both of these two states. Specifically, the ra- 
tio of the the intercavity coupling and on-site interaction 
strength for the transition from Mott-like region to SF- 
like region is modified by a time dependent increment. 
Moreover, the superfluidity of a initial SF state will break 
owing to the decay of the off-diagonal long-range order. 

Consider a system consisted of atoms and cavities cou- 
pled weakly to a bosonic environment at zero tempera- 
ture. As the dimension of individual cavities is generally 
much smaller than their spacing, we assume the photons 
emitted from each cavity are uncorrelated. The total 
Hamiltonian therefore reads 



H — Hs + Hba- 



th 
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k '^k^f^j fi,ka, the Hamiltonian for environment. 



and Hcoup = T,iT,a,kiv*k 
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aj + h.c.) the coupled 



•-coup — /^j /^a,k\'lka' j,ka^3 

term, a — a,c labels the operators and physical quanti- 
ties associated with atoms and cavities, respectively. iOk^ 



denotes the frequency of environmental modes. 



and 



Tj^fc^the creation and annihilation operators of quanta in 
the kath model at site j, and rjk^ the coupling strength. 
Here we set H = 1. 

The system we modelled, as depicted in Fig. 1, 
is a two-dimensional array of resonant optical cavities, 
each embedded with a two-level (artificial) atom coupled 
strongly to the cavity field. The possible realizations 
such as photonic bandgap cavities and superconduct- 
ing stripline resonators et al. With and uJc be- 




where Hs is the Hamiltonian for the system, Hbath 



bare dressed 
basis basis 



FIG. 1: A type of possible topologies for two-dimensional 
cavity arrays, (a) Individual cavities are coupled resonantly 
to each other due to the overlap of the evanescent fields. Each 
cavity contains a two-level system coupled strongly to the cav- 
ity field and immerses in a bosonic bath(marked by the dash 
line), (b) Energy eigenvalues of individual cavity-atom sys- 
tem on each site, cjc ~ ijJa is assumed for simplify. The anhar- 
monicity of the Jaynes-Cummings energy levels can effectively 
provide an on-site repulsion U to block the absorption for the 
next photon. 
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ing the frequency of atom transition and cavity mode 
respectively, in the rotating wave approximation(RWA), 
such individual atom-cavity system on site j is well de- 
scribed by the Jaynes-Cummings Hamiltonian, H^'-^ = 

(jJaOjCij + LUcCjCj + PiojCj + h.c ). Here a| and aj(c|, Cj) 
are atomic (photonic) rasing and lowering operators, re- 
spectively, (3 the coupled strength. In the grand canoni- 
cal ensemble, Hs is therefore given by combing ^//"^ with 
photonic hopping term and chemical potential term, 

^« = E - E ""h' - E '^"j- ■ (2) 

Kjji is the photonic hopping rate between cavities. Since 
the evanescent coupling between cavities decreases with 
the distance exponentially, we restrict the summation 
j') running over the nearest-neighbors, rij — a^^aj + 

cJ^Cj is the total number of atomic and photonic excita- 
tions on site j. /i is the chemical potential, where the 
assumption /i — fij for all sites has been made. 

Due to the strongly coupling, as shown in Fig. 1(b), 
the resonant frequencies of individual atom-cavity system 
are split into 



E\±.n) = nuJc± yn/S'^ + — ~ —, (3) 

where |±,n) labels the positive(negative) branch of 
dressed states, A = uJc — uja is the detuning. The an- 
harmonicity of the Jaynes-Cummings energy levels can 
effectively provide a on-site repulsion. For instance, the 
resonant excitation by a photon with frequency ^'|±,i) 
will prevent the absorption of a second photon at i\, 
which is the striking effect known as photon blockade ^ . 
It is therefore feasible to realize a quantum simulator in 
terms of the system described by Eq. (2). This so called 
Jaynes-Cummings-HubbardfJCH) mode is recently sug- 
gested by Greentree et al. [2]. 

However, the situation changes dramatically once tak- 
ing the coordinates of environment into consideration, as 
described by Hamiltonian(l). A non-equilibrium dynam- 
ics for open quantum many-body system do arise, which 
is a formidable task to solve. Here we propose a treat- 
ment to eliminate those external degrees of freedom. To 
approach this, we rewrite Hamiltonian(l) as 

H = Hiocal - E '^Jf4^r - E /^'^J ' 

where ff/oca/ = J2j Hf'^ + Hbath + H^oup- 

First considering the case that the jth cavity contained 
a initial photon interacts with a bath, the dynamics is 
governed by 

Hj ^ ^^cc]cj +^^kj]^kj3.k. + ^i^l/j^k^Cj + h.c). 
k k 

(5) 

We denote its eigenvalue as w and expand the eigenvec- 
tor as = ecc]|0) + J2k'^krlj9). Cc and Cfe are 



the probability amplitudes for the excitation occupied by 
cavity field and environment, respectively. |0) denotes 
the vacuum state. Deducing the equations of these two 
amplitudes, one can express ek in terms of Cc and, un- 
der the Born-Markov approximation, integrate out the 
degrees of freedom of environment, and obtain 

(wc + 5ujc - i^c)ec = iocc- (6) 

Sujc is known as an analog to the Lamb shift in atomic 
physics and significantly small when the coupling to en- 
vironment is weak. 7c is the decay rate and indicates a 
finite lifetime of cavity mode [l5| . 

This motivates us to introduce a quasi-boson described 
by Cj with a complex eigenfrequency fie = i^c—ijc-, where 
Sujc has been absorbed into oJc, to effectively describe the 
open system described by Hamiltonian (5) in terms of 

i?f^|</),)=c.f^|^,), (7) 

with the effective Hamihonian H^-^-^ — QcClCj and now 

— ecCj|0). Because of loss, the system would be 
nonconservative and corresponding operators would be 
non-Hermitian. The communication relation of Cj reads 
[Cj,Cj,] = (1 + i^)Sjj'. Recognizing ^ is in order of 
^, with Q being the quality factor of individual cavity. 
The bosonic communication relation is therefore approxi- 
mately satisfied for the high-Q cavity, which is meeting in 
most experiments about cavity quantum electrodynam- 
ics(QED). 

The complex eigenfrequency underlines the facts that, 
on one hand, dissipation is the inherent property for re- 
alistic cavity. When a photon with certain frequency has 
been injected into a dissipative cavity, the composite sys- 
tem of cavity-filed plus environment cannot be charac- 
terized merely by the frequency of the injected photon, 
however, we must take the impacts of environment into 
account. On the other hand, in general we are not con- 
cerned the time evolution of bath. In this way, the array 
of dissipative cavities can be regarded as a configuration 
consisted of quasi-bosons. Quite similar operations can 
be performed on atom to introduce another kind of quasi- 
boson described by Aj with the frequency fla = lo^ — ija, 
where 7a is the atomic decay rate. 

We can therefore rephrase Hamiltonian (1) with the 
renormalized terms, 

^ = E<^- E ^n'C}Q,-Y,fm,, (8) 

J <jj'> 3 

with now Hf^ = naA]Aj+ncC]Cj+/3{A]C.j+h.c.) and 

rij = A]^Aj + CjCj. One key feature of Hamiltonian (8) is 
now the losses describe by leaky rates 7a and 7c but not 
by field oscillations. Without having to mention the ex- 
ternal degrees of freedom, this effective treatment would 
be of great conceptual and, moreover, computational ad- 
vantage rather than the general treatment as Hamilto- 
nian (1). A more microcosmic consideration points out 
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that, in cavity QED region, since the atom is dressed by 
cavity field, the atom and field act as a whole subject to 
a total decay rate F [l^. In particular, F = n(7a + 7c) 
for A = 0. 

To gain insight over the role of dissipation in the SF- 
Mott phase transition, we use a mean field approximation 
which could give reliable results if the system is at least 
two-dimensional We introduce a superfluid param- 
eter, -0 = Re{Bj) = Re{Bj). In the present case, the 
expected value of Bj {Bj ) is in general complex with the 

formation Bj = ip — iip^{B'lj = ip + iip^). is a solvable 
small quantity as a function of 7q and 7c, and vanishes 
in the limit of no loss. Using the decoupling approxima- 
tion, b]Bj, = {B])By + {Bj,)b] - {B]){b]), the result- 
ing mean-field Hamiltonian can be written as a sum over 
single sites, 

jjMF ^ Y,{H;"-z^m]+B,)+z^m'-^in,+o{^l)}, 

(9) 

where we have set the intercavity hopping rate Hjji — 
K for all nearest-neighbors with z labelling the number. 
Corresponding to the order parameter in related closed 
system, whether -0 vanishes or has a finite value can be 
used to identify the Mott-like and the SF-like regions. 

ip can be examined analytically in terms of second- 
order perturbation theory, with respect to the dressed 
basis. For energetically favorable we assume each site is 
prepared in the negative branch of dressed state. But be- 
cause the dressed basis is defined on n > 1, a ground state 
|0) with the energy i?|o) = need to be supplemented. 
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(10) 



X and Q arc functions of all the parameters of the whole 
system. Since the evanescent parameter k is a typical 
small quantity in systems of coupled cavities, the per- 
turbation theory gives good qualitative and quantitative 
descriptions comparing to the numerically results given 
by explicitly diagonalizing 0, Ull . 

Arguably the most interesting situation is the photon- 
photon interactions are maximized, namely, cavities on 
resonant with atoms and with one initial excitations per 
site [li|. In addition, T = -fa + Ic = 1- With Fi = 
ujc — (3 — ^J, and F2 — —ujc + (\/2 — 1)P + /i, in eq. (10), 



e 



3+272 ^ n J 

4Fff4^ > and 



F, 



X 



(3 + 2V^)F2 1 



2F^ + 272 4F2 + 472 



-2jt ■ 



(11) 



In the absence of loss, one can recognize x = is the 
well known self-consistent equation and therefore distin- 
guish the SF phase and Mott phase. Nevertheless, the 
coupling to environment inducing a non-equilibrium dy- 
namics, and there is thus no strict phase exist. How- 
ever, there remains a time-dependent boundary given by 
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FIG. 2: The transition from Mott-like to SF-like region. In- 
fluences of dissipation on the critical ratio depend on the leaky 
rate 7 (the dot and solid line for ^ = 0, 0.05, respectively) 
and the time one reaches the transition (the red, green, and 
blue lines for t — 0, 0.l7~^, 0.27~^, respectively). 



-0 which can distinguish two slowly decayed and quali- 
tatively different regions, i.e. the SF-like and Mott-like 
regions. 

In what follows, we discuss the physics of the transition 
between these regions from two aspects. First, we start in 
the Mott phase and discuss the impacts of dissipation on 
transition ratio (^)c- Consider the initial state is deep in 
the Mott phase, ^ — 0, and we continually increase the 
intercavity coupled rate. For the related ideal case, we 
wiU reach the SF phase at ^ = (^)c ~ 0.16. However, 
in the presence of dissipation, owing to the continuous 
leakage of coherence the effective tunnelling energy will 
be lower than what we set. Stated in another way, the 
photon hopping rate not only dependent on the coupling 
parameter between cavities, but also on the coherence in- 
side cavities. As shown in Fig. 2, one cannot expect the 
appearance of photon hopping at ^ ~ 0.16. We must 
continuously increase k and the rate of increase must 
faster than the decay of coherence. Only in such way can 
the transition occur before the system reach the steady 
state. However, even we start exactly at ^ = (■^)c, 
we also cannot expect we are in the SF state. As the 
dissipation is the inherent nature of open quantum sys- 
tem, the critical ration (^)c, even at t = 0, is modified 

by an increment ~ -p-. This can be understood in a 
simplified picture. If there was no excitation in system, 
i.e. we are in |0), the decay would not happen. However, 
the open quantum system is inherently different from the 
correspondingly perfected system. 

In contrast, we start with the SF phase in the vicinity 
of transition ratio and track the time evolution of sys- 
tem. The exponential term in Eq.(lO), e"'"*, indicates 
the expected decay of ij) and related physics quantities. 
However, the more important impacts of the presence of 
external environment are revealed by x- As illustrated 
in Fig. 3, for t ^ ip has a slightly decrease scaled 
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FIG. 3: The temporal decrease of the superfluidity and the 
photon number fluctuation on each site for a certain initial 
state(inset). For a initial SF state (the red line for ^ = 0.3 
and blue line for — 0.2), before tc the long-range order 
decays continuously and the fluctuation on each site is a to- 
tal effect of photon hopping and photon leakage (the solid 
line for j = 0.01 and dash line for j = 0.02). Beyond tc, 
the superfluidity breaks down and the related photon num- 
ber fluctuation behaves as the fluctuation of a Mott-like state 
(dot-dash line). 



by ?2- For t 



> (3 ^, the leader term is zne It is 



this term pronounces the decrease of effective tunnelling 
energy which, mathematically, originates from the non- 
diagonal element in Hamiltonian (9) respecting to the oc- 
cupation number basis. Consequently, a photon hopping 



rate given initially in SF phase will cross the critical 
point at a time tc — Before t^, a non-local re- 

gion is still recognized as non-local. The dissipation does 
not change the fundamental nature of the system, albeit 
with the decrease of long-range coherence. Nevertheless, 
beyond t^, the superfluidity breaks down and gives rise 
to a localized wavefunction, i.e. a transition to Mott-like 
region do occur. An analogous effect is recently described 
in optical lattice system, where a time-dependent tunnel- 
ing rate controlled by external filed leads to a sweep from 
the SF to the Mott phase [13]. Differently, in our case 
the decay of long-range order is resulted from the leak- 
age of coherence, and the total number of excitations is 
not conservative. As a result, there remains a statistical 
fluctuation acted on each site. 

In conclusion, we have shown analytically the features 
of SF-Mott phase transition in the coupled cavity arrays 
in the presence of dissipation. Our analysis fully takes 
into account the intrinsically dissipative nature of open 
quantum many-body system, and identifies how dissipa- 
tion and decoherence would come into play. For further 
experiment realization, we predict there is a localized ten- 
dency characterized by the suppression to the superfiu- 
idity. 
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